Melting and Dimensionality of the Vortex Lattice in Underdoped YBa 2 Cu 3 Oe 
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Muon spin rotation (/iSR) measurements of the magnetic 
field distribution in the vortex state of the oxygen deficient 
high-r c superconductor YBa2Cu306.eo reveal a vortex-lattice 
melting transition at much lower temperature than that in the 
fully oxygenated material. The transition is best described 
by a model in which adjacent layers of "pancake" vortices 
decouple in the liquid phase. Evidence is also found for a 
pinning-induced crossover from a solid 3D to quasi- 2D vor- 
tex lattice, similar to that observed in the highly anisotropic 
superconductor Bi2+ :c Sr2-a;CaCu2 08+2 / . 
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In the high-temperature superconductors there exist 
several exotic vortex lattice (VL) phases owing to the 
weak coupling between the superconducting Cu0 2 layers 
(which gives rise to highly flexible vortices), a relatively 
short coherence length (which enhances their susceptibil- 
ity to pinning) and high values of T c (which permit large 
thermal energies to be reached in the vortex state). To 
date much attention has been paid to VL melting and 
pinning effects in the highly anisotropic superconductor 
Bi 2+2 ,Sr 2 _ ;c CaCu 2 8+ < / (BSCCO) and in YBa 2 Cu 3 7 -,5 
(YBCO) for near optimal oxygen concentrations (i.e. 
<5«0.05, giving the highest value of T c ). The vortices in 
BSCCO are often described as stacks of two-dimensional 
(2D) "pancakes", which can become misaligned at ele- 
vated temperatures and/or in the presence of strong ran- 
dom pinning sites. On the other hand, in "optimally" 
doped YBCO the vortices behave in a three-dimensional 
(3D) manner. This difference is attributed to the mass 
anisotropy 'y—(m c /m a i i ) 1 ' 2 ~X c /X a i ) Ki5-7 in YBCO and 
7 « 50-250 in BSCCO jy, where A a b and A c are the mag- 
netic penetration depths describing the screening of flux 
by supercurrents flowing in and out of the Cu0 2 layers, 
respectively. An added complication in many studies of 
YBCO is the presence of twin-plane boundaries which 
may act as extended pinning sites for vortices. 

The muon spin rotation (/iSR) technique measures the 
muon precession frequency distribution and thus the in- 
ternal magnetic field distribution n(B), also known as the 
/iSR line shape, which has proven to be a powerful tool 
for investigating VL phases 0-0]. In optimally doped 
YBCO, strong local pinning broadens the /iSR line shape 



II as predicted for small random displacements of 3D 
vortex lines p|. By contrast, random pinning in BSCCO 
leads to a field-induced reduction of the /xSR line width 
|g,0-||] associated with a dimensional crossover from a 
3D-VL to a quasi-2D system consisting of independent 
weakly interacting VLs in different Cu0 2 layers. In the 
quasi-2D region, the repulsive interaction between pan- 
cake vortices in the layers exceeds the strength of their 
interlayer coupling, so that random pinning results in a 
misalignment of the pancake vortices in the field direc- 
tion. 

The effect on the /iSR line shape of thermal fluctua- 
tions of the vortex positions has been the focus of nu- 
merous studies in BSCCO ]2|, pi]| , pd| , [l3| . As explained in 
Rcf. PCI , the typical time scale for vortex fluctuations 
is short enough that the muon detects a fluctuation- 
averaged field. This results in a reduction of both the 
width and the skewness of the /iSR line shape. In clean 
samples the VL melts in a first-order phase transition 
(see, for example, Ref. Jl2|]). Recent /xSR measurements 
in BSCCO have been interpreted as evidence for a two- 
stage VL transition: first the intralayer coupling of the 
pancake vortices is overcome by thermal fluctuations, 
then their interlayer coupling is lost |13| ] . 

In the less anisotropic compound YBCO, a melting 
transition has been observed by small-angle neutron scat- 
tering, magnetization, transport, specific heat and Hall 
probe ac susceptibility measurements (see, for example, 
Ref |14|). However, unlike in BSCCO, the vortex- liquid 
phase is found only in a very narrow region below the 
second-order phase transition at B c2 (T). A melting tran- 
sition near B c2 (T) is difficult to study with the /xSR tech- 
nique, because even for a 3D VL, the /iSR line shape nar- 
rows and becomes more symmetric due to the overlap of 
the vortices. 

It is well established that removal of oxygen from 
the CuO chain layers in YBCO weakens the coupling 
between the superconducting Cu0 2 layers. Here we 
report /iSR measurements of the VL in the oxygen- 
deficient compound YBa 2 Cu3O6.60; which has an in- 
creased anisotropy ratio 7 « 22-36 (see, Ref. [|15| ) . We 
observe an expanded vortex-liquid region and a pinning- 
induced crossover from a 3D to a quasi-2D system, quali- 
tatively resembling the vortex phase diagram of BSCCO. 

The fiSK experiments were performed on the M15 and 



M20 surface muon beamlines at TRIUMF, Canada, us- 
ing the experimental setup described in Ref. Pq| . Mea- 
surements of the internal magnetic field distribution were 
made in both twinned and detwinned crystals of under- 
doped YBa2Cu3O6.60i with superconducting transition 
temperatures T c — 59(0.1) K. These samples were pre- 
viously studied with /iSK at low temperatures Q]. The 
crystals were mounted with their c-axis parallel to both 
the applied magnetic field Hand the muon beam direc- 
tion. The positive muons were injected into the sample 
with their initial spin polarization perpendicular to H . 
As described fully elsewhere JlT]], n(B) or the muon pre- 
cession frequency distribution n(u> = 7^-B), where 7^ is 
the muon gyromagnetic ratio, is obtained by monitoring 
the e + count rate as the muon decay pattern sweeps by 
the positron detectors. 

Figure 1 shows the evolution of the fast Fourier trans- 
form (FFT) of the muon precession signal upon warming 
the sample, following field cooling to T = 2.4 K in a mag- 
netic field #=1.49 T. At T = 20.5 K the FFT shows the 
basic features expected for a well-ordered 3D solid VL 
in particular a high-field "tail" associated with the 
region close to the vortex cores. However, at T = 40 K 
the FFT is completely symmetric and the line width is 
drastically reduced -- both of which characteristize a 
vortex-liquid phase. Despite these obvious changes in the 
/iSR line shape, the melting transition cannot usually be 
determined accurately by visual inspection, because the 
output frequency spectrum is artificially broadened by 
the finite time range and the "apodization" needed to 
eliminate "ringing" in the FFT |18|. To quantify these 
changes in the field distribution we calculate the skew- 
ness parameter H, a = ((SB) 3 ) 1 / 3 / ((SB) 2 ) 1 / 2 , where 
((SB) n ) = ((B-(B)) n ). To obtain reliable values for 
the moments ((5B) n ) we fit the nSR spectra in the time 
domain to a polarization function calculated assuming 
a Ginzburg-Landau (GL) model for the VL field profile 
B(r), as described in Ref. M. The fitted function B(r) 
properly accounts for the long high-field "tail" in n(B) 
when the VL is 3D. This method of analysis also provides 
a means of monitoring the /iSR line width through the fit- 
ted value of 1/A^ b . A hexagonal arrangement of vortices 
was assumed, consistent with recent neutron scattering 
measurements on untwinned YBCO p9|. 

Figure 2 shows the temperature dependence of a and 
1/A^ b at three of the fields considered. The observed 
values of aw 1.2 are in agreement with values previously 
obtained in the 3D-VL phase of BSCCO [§|]|{Q]I! and 
that predicted for s- and ci-wave VLs at low reduced fields 
B/B C 2 [£0|. For temperatures above T m the fitted val- 
ues of 1/X 2 b and hence a drop to zero. This does not 
mean that the superfluid density has decreased to zero, 
but rather that the £tSR spectrum no longer contains a 
component corresponding to an ordered 3D vortex struc- 
ture. For T > T m the /iSR spectra fit well to a single 
Gaussian function, with a line width much greater than 
that due to nuclear dipolcs. 

In Ref. |gl| a sharp transition was observed in ac 



susceptibility measurements on twinned YBa2Cu3Og.60 
crystals below a characteristic field H* «0.07 T, possibly 
associated with a melting transition. Above H*, where 
the transition was found to be continuous, it was sug- 
gested that a quasi-2D system exists which is highly sen- 
sitive to pinning-induced disorder. In the present study, 
where H > H*, the field distributions for T < T m and 
0.1 < If < 1.49 T fit well to the GL model for an ordered 
3D VL. The rms deviation of the vortices from their ideal 
positions in the hexagonal VL was found j| to be < 8 %. 
These results suggest that the VL for our samples in this 
region of the B-T phase diagram exhibits a 3D behav- 
ior, most likely consisting of stacks of strongly coupled 
pancake vortices. Above T m , thermal fluctuations of the 
vortices lead to a loss of long-range spatial order. 

The top panel of Figure 3 shows the //SR line shape in 
the twinned sample of YBa2Cu3Og.60 after cooling in a 
field ff«2.89 T to T = 2.5 K, followed by an increase of 
AH w 0.01 T. Notice that the small background signal 
is positioned at the external field H = 2.90 T, whereas 
the signal originating from the sample looks as though 
the external field is still H = 2.89 T. This implies that 
the VL is strongly pinned. However, unlike in the line 
shapes observed at the lower fields, the "tail" appears 
on the tow-field side of the cusp. Numerical calculations 
performed in Ref. [J22fl , which account for the sample ge- 
ometry effect, show that such a line shape can originate 
from a system of 2D pancake vortices that are ordered 
in the planes but uncorrelated between adjacent layers. 
The low-field "tail" is specifically associated with a lower 
flux density at the sample edges due to a nonuniform 
demagnetization. Thus clear evidence is observed for a 
pinning-induced dimensional crossover above a critical 
field B a (T). 

By contrast, at the same field H = 2.9 T in the de- 
twinned sample the "tail" extends to the high-Held side 
of the distribution, as expected for an ordered 3D vor- 
tex structure. It is tempting to attribute the opposite 
skewness of the /iSR line shapes in the two samples to 
the presence of twin planes. However, because the twin 
planes extend the full depth of the sample, they would 
displace the 2D VLs by the same amount in all layers. 
Since the two samples were not from the same growth 
batch, the difference is more likely related to the concen- 
tration and randomness of pointlike defects. This result 
stresses the sensitivity of the VL structure to disorder at 
this field. 

The middle panel of Fig. 3 shows that the vortices 
begin to depin at T p « 19.8 K due to thermal fluctu- 
ations, and reposition themselves with an average field 
5 « 2.90 T at T = 23 K (bottom panel of Fig. 3). Addi- 
tional fiSK spectra taken without shifting the field show 
that the VL melts at T m «20 K. 

Figure 4 summarizes our results in a phase diagram. 
The upper critical field line represents an approximation 
assuming B c2 (T)=B c (0)[l-(T/T c ) 2 } and £ c2 (0) = 70 T. 
The data for the twinned and detwinned samples appear 
to fall on the same melting line. A fit of the data below 



if = 2.9 T to the phenomenological melting curve [O 
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yields to = 1.6(1), where X ab = A o6 (0)[l - (T/T,) 4 ]- 1 / 2 , 
and K and to are constants. This is well below either 
the predicted value to = 4 when interlayer coupling of 
the pancake vortices is dominated by electromagnetic in- 
teractions, or to = 3 when Josephson coupling becomes 
important p3[. As discussed in Refs. fOJ23fl, for finite 
Josephson coupling an exponent to = 2 is expected for an 
additional thermodynamic transition at B < B cr which 
decouples the layers in the liquid phase. The decoupling 
line is predicted to be p3] 
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where Cd ~ 0.1 is a decoupling constant and s is the 
interlayer spacing. The dashed curve in Fig. 4 is a fit 
to Eq. (|). Taking s = 11.7 A ||], A Qb (0) = 1642 A |] 
and 7 = 22-36 Jl|], we calculate c D = 0.025-0.067 from 
the fit. A value Cd =0.076 was obtained for BSCCO in 
Rcf. |ll|] . To determine whether the melting transition 
(i.e. to a liquid of vortex lines) coincides with the low- 
field interlayer decoupling transition, more measurements 
in the vicinity of the phase transition are needed. We 
note that, contrary to the conclusion in Ref. p[ , recent 
detailed /^SR measurements in BSCCO fl3| suggest that 
B m (T)<Bl(T). 

In conclusion, we have observed changes in the fiSK 
line shape of underdoped YBa2Cu3O6.60 which identify 
a VL melting transition far below that of optimally doped 
YBCO. More precisely, fits to the data suggest that there 
is a decoupling transition from a liquid of vortex lines to 
a liquid of 2D vortices. Our measurements also establish 
the existence of a pinning-induced dimcnsionl crossover 
to a quasi-2D vortex system, similar to that observed in 
the highly anisotropic material BSCCO. 
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FIGURE CAPTIONS 

Figure 1 . Fourier transform of the muon precession sig- 
nal in twinned YBa 2 Cu 3 O 6 .60 at T = 20.5 K and T = 40 K, 
in a field H= 1.49 T. 

Figure 2. Temperature dependence of the skewness a 
[top] and 1/A 2 fc [bottom] in detwinned YBa 2 Cu3O6.60 at 
H = 1.25 T [crosses] and in twinned YBa2Cu3O6.60 at ap- 
plied magnetic fields H = 0.74 [open circles] and 1.49 T 
[solid circles]. 

Figure 3. Fourier transforms of muon precession sig- 
nals in twinned YBa2Cu3O6.60- Top panel: after field 
cooling atJTw 2.89 T to T = 2.5 K followed by an in- 
crease in the field to 2.90 T. Middle and bottom panels: 
after subsequently warming the sample to T= 19.8 K and 
23.0 K, respectively. 

Figure 4. The vortex B-T phase diagram for 
YBa 2 Cu3O6.60- Open and solid circles correspond to 
the twinnned and detwinned samples, respectively. The 
dashed curve B^ C (T) is a fit to Eq. (g), the solid curve 
B C 2 (T) is a theoretical line for the upper critical field (see 
text) and the solid curve B m (T) is a fit of all of the data 
to Eq. (|l|). The transition line B cr (T) is approximate. 
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